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Two one-dimensional Luttinger liquid systems coupled by Hund's coupling are studied by 
the renormalization group and the non-abelian bosonization methods. It is found that the Hund's 
coupling is always relevant irrespective of the repulsive interaction between electrons within each 
chain. The properties of the resulting strong coupling fixed point are discussed. 
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I. INTRODUCTION 

One-dimensional (ID) electron systems provide examples where the electron correlation effects play a dominant 
role. The metallic phase of ID interacting electron system is characterized by the separated gapless charge and spin 
excitations, and it is called the Luttinger liquid phase. Recently there has been great interest on the two coupled Lut- 
tingeiiiijquid systems as the model for ladder systems, zigzag system, stripe phase in cuprate superconductor, and nano 
tubesBH. Theoretically, the relevancy of inter-chain hopping t±U, the opening of spin gap, and the superconducting 
instabilities have been studied.QT] 

In this paper, we investigate two ID Luttinger liquid systems coupled by [Hund's coupling without interchain hopping 
t±. This coupled system is interesting as a first step towards 2D Sr2Ru04l3'tl. AiHubbard chain away from half-filling 
can be modelled as a Luttinger liquid for arbitary values of on-site repulsion f/Ei, and two coupled Hubbard chains in 
the weak coupling regime U, Jr < t can be investigated using the bosonization technique {t is a band width, Jh is 
Hund's coupling). Nagaosa and OshikawaQ investigated two coupled Luttinger liquids using a semiclassical analysis 
under the condition 6 max{t,t±) < \Jh\, where 6 is the density measured from half- filling. They showed that a spin 
gap opens up in this system and that the charge fluctuations between chains are massive. For generic values of 
t, t± , Jh , their result is valid only near half-filling. 

In present paper we investigate this model at arbitrary filling factors away from half filling. We focus on the influence 
of Hund's coupling on the charge and spin excitations and its interplay with the repulsive Hubbard interaction. In 
order to preserve SU(2) invariance in renormalization group analysis we use non-Abelian bosonization. The weak 
coupling phase diagram can be mapped out from the analysis of the resulting perturbative R.G. equations. We 
find that Hund's coupling is relevant even for infinitesimally small initial value of Jh, and the system fiows into a 
new strong coupling fixed point. At the strong coupling fixed point, all spin excitations are gapped. As for the 
charge exciations, the symmetric charge excitations 4'p+, 9p+ (See Eq.(|9|))remain gapless, while the asymmetric charge 
excitations (pp^, Op- are gapped. Thus, our result extends the results by Nagaosa and Oshikawa to general filling for 
the case t± — 0. Furthermore, we have investigated the dependence of the physical properties on the value of the 
Hubbard repulsion, which is not addressed by Nagaosa and Oshikawa: We find that, even the system is in gapped 
phase for all < < t and < Jh < t, the gap (i.e. the scale where the renormalized coupling constants diverge 
or become the order of the bare cut-off scale) strongly depends on the relative magnitude of U, Jh ■ U > Jh , the 
strong coupling fixed point can be reached faster than the opposite case Jh > U. 

This paper is organized as follows: In section II, we specify the model and re-cast it in a bosonizcd form. In section 
HI, the renormalization group equations are derived using the operator product expansion method. We present the 
renormalization group flows in section IV and conclude with discussions in section V. 

II. MODEL 

We consider two Hubbard chains coupled by Hund's coupling. 

H = -t^cj^^cija + U^Ui^nui ~ ^ JhSu ■ S2i ^ Hq + Hu + Hh, (1) 

i,l i 
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where a is a spin index and Z = 1, 2 is a band index. Jh is positive. Note that there is no inter-chain hopping. We 
assume the imphcit normal ordering of all interaction terms in the Hamiltonian. For simplicity, the band dependences 
of hopping parameter tij^i and Hubbard interaction Ui are ignored, two coupled chains are considered instead of three, 
and only the nearest neighbor hopping is considered. 

We investigate only the weak coupling regime U < t, Jh < t, in which the linearization of non-interacting electron 
spectrum near Fermi points and the pcrturbative treatment of interactions are legitimate. A generic filling away 
from half-filling is considered, so that the umklapp process can be neglected. Because the inter-chain hopping term is 
absent, the non-interacting energy spectrums of two bands are degenerate. The linearization gives the decomposition 
Ciail \fa ~ TpRaiix) e*'^^^ -I- tl^Lai{x) e"**-'^^. a is a lattice spacing. Substituting the decomposition into Hq of Eq.(P, 
we get 

-ff = ^ / dxVF (V'ffla ^dxIpRla " TpUa ida^TpLla^ , (2) 

where vp = 2tasm{kFa)- To express the Hamiltonian in the bosonized form, it is convenient to introduce the (chiral) 
charge and spin current operators. 

jR,l^^i}R^l^ai^R,l,a, ^R.l=^^R,l,a^^R,l,l3- (3) 

a af3 

The left masking currents are defined analogously. The operator product expansion allows us to express Hq in terms 
of currentsEJ. 

Ho = Y.j dx\^[Jh + Jh) + ^(Jl; • iLi+iRi ■ im)] ■ (4) 

The Hubbard part Hjj can be expressed in terms of currents alone, while the Hund's coupling term Hh cannot be, 
since 

S,/a = J,, + J;,, + e-2^'=- i^m + Ai.^ i^RW- (5) 

At this point, we apply Abelian and Non-abelian bosonization to the charge and—spin part, respectively. The phase 
field (jipi and its conjugate momentum Hp; are introduced for the charge currentsM. 

Jr.i + Jli = ]J^dx4>ph Jri - Jli = -\j-^^pi- (6) 
And for the spin currents SU(2) matrix field gi{x, r) is introduc© 

J£; = --tr(a,5;gr'v)' Jri = +-^^{97'd-z9i^)^ (7) 

with z — X -\- ivpT, z — X — ivpT. The elements of the matrix field gi can also be expressed in terms of scalar field 
through the vertex operator constructionlU 

(t)„i is the phase field of the Abelian bosonized spin currents, and O^ri is the conjugate field. The above Abelian 
bosonized representation is useful in interpreting the results of R.G. equations physically. For later convenience, we 
define the symmetric and asymmetric charge and spin modes. 



4>p± = ^[(t>pl±(l)p2j, 0<T± ± (/)ff2 j ■ (9) 

According to the Non-abelian bosonization rulc0"lii the fcrmion bilinear can be represented as 

T.^^^^^^^ - tr(5/(-)f ) e-^^-(-). (10) 

a/3 



2 



As a result of the parametrizations of currents and fermion bilinear Eq.(|6|,|7|,|10D, it is possible to express the total 
Hamiltonian in terms of the scalar fields (f>pi and the SU(2) matrix fields gi. For the renormalization group analysis, 
the Euclidean Lagrangian formulation is more convenient. 



d6 



•pi 



dVcT 



dx 



wzw 



Ai E / (IxdrSm ■ J 
Ixdr \\2 ( 



Ll 



- / da;dr A2(Jli • Jl2 + J_ri • Jm) + ^zi^Li ■ J_R2 + Jfli • J 



L2) 



A4 / dxdr 



(11) 



As / dxdr 



tr(5l(a;)|) e'^^^-i^") tr(.g2(a:)|) e+'^'^-^(") + h.c 
[tr(5j(a;)) e^*^"*''!^^) tr(.92(a:)) e+'^'^''2(^) + h.c 



2/^2 



E / dxdT2 



'pi d(j)p2 \ ( d(j)pi dcj)p2 



dvcT dvr 



dx dx 



(12) 



(13) 



where 



and 



VF \ I 



Ua 

TTVp 



(14) 



if. = iA/i + ^. 

TTVp 



(15) 



wzw 



is the so-called Wess-Zumino-Witten action, and we don't need its explicit form. We just note that the 



Fermi velocity for the matrix field g is renormalized by the interaction vp Vg 
(bare value) of Ai is 2Ua > 0, and the initial value of A2, A3 and A4 is Jhci > 0. 



Vpil 



Ua 



). The initial value 



Although the initial values of A2, 
A3 and A4 are identical they will scale differently under the R.G. flow. The I/K2 and A5 terms in the above action is 
absent in the bare action (the initial values are 0), but they are generated by the second order perturbation. 

The action Eq.(p|) defines a critical fixed point of two independent U{1) Gaussian model and two independent 
SU(2) fc = 1 Wess-Zumino-Witten (WZW) theory. The terms in Eq.(|l|)(Ai, i 1, 2, 3, 4, 5) are the perturbations to 
the fixed point. We study how the perturbations influence the fixed point defined by Eq.(pT[) via the renormalization 
group method in the next section. 



III. RENORMALIZATION GROUP EQUATIONS 



It is well-known that the renormalization group (R.G.) equations for the fixed point Hamiltonian perturbed by a 
number of scaling operators, up to the second order in couplings, are determined by the scaling, dimensions of the 
scaling operators and the operator product expansion (OPE) coefficents of the scaling operatorstj. Explicitly, 



dXi 



d \nL 



= (d- Xi)Xi - E CijkXjXk H : 

jk 



(16) 



where Ai is the coupling constants of the scaling operator Oi,d=l + l = 2is the space-time dimension, and xi is the 
scaling dimension of the scaling operator. L is a cut-off length or time scale. The OPE coefficents c^fc is determined 

by 



o.{r)oM-T.-Sh-om 



(17) 



Among the perturbations in Eq.(12), only A4 has a scaling cliH iengip n different from 2, which is a;4 = 1 + Kc 
The relevant operator product expansions for our case arell3^t£HlZl 
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gia(t> pi {x,T)^-ia(t> pi (0,0) 



1 + za(z,5,>p, (0, 0) + <^p/ (0, 0)) 



with Zc = X ^ iVcT, Zc 



X — IVrT. 



The OPE of WZW model (fc = 1) are given by 

1 k ShnSab 



Jli{z){x,T)g„,{Q,Q) 



1 ie'" 



be 



Olrn 

27rz, 



tr 



{giix, T)a'^) tr(gtjO, 0)^") ^ \Sab + *e'^''^(z. J£,(0, 0) + J]^,(0, 0)) 



+ ■ 



(18) 



(19) 
(20) 

(21) 



where Zg = x + iVsT, Zg — x — iVsT, and t° are the generators of SU(2) algebra in the spin basis t°- = ti{t°-t^) = 

[f", t''] = it"-^'^ f^. The terms of dimension higher than 2 are omitted since they are irrelevant in our computation. 
The OPE Eqs.(|l9[|2^) are the fundamental consequencs of the chical gauge invariance of WZW model. The OPE 
Eq.(pTl) follows from the fusion rule of Kac-Moody primary fieldfilil. For level 1 (fc = 1) S1I(2) WZW theory, the 
product of spin 1/2 primary fields can only give rise to the identity field and its descedantslljo. J and terms in 
Eq.(^I|) are the descendant fields of the identity field. 

Using Eqs.(^6|jl8|jl||^,|l|), the R.G. equations can be computed straighforwardly. The results are : 



dO 



d\nL TTVcVs^^^Vs^' 7o (cos2 + u2sin2 0)^<=-i' 



(22) 



d\nL 



nVcVs Vs 



(23) 



dXi 

d\nL 



d\2 



2L 

2t:Vs 



d\nL 



= 0. 



dXs 
dlnL 



27r V V, 



(24) 
(25) 
(26) 



dXi \ \ , ^l'^4 -^3^4 

= (1 - Kc)M 



d\nL 
dX5 



AtTVs 



(1 - K,)X, - 



2'KVs 



■ A3 A4 . 



(27) 



(28) 



dlni IGnVg 

We note the similarity of the above R.G. equations with those of ID Kondo lattice problem0. The Eq.( |2^j2^ ) suggests 
a combination of fields (j)p± — -^{4' pi i 0p2)- In terms of (j)p±, the sum of the charge parts of our action becomes 



2 ^Kc K2' 



dxdr 
dxdr 



( dVcT 
\ dVcT 



dx 

d(l>p- 
dx 



(29) 
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Then from Eq.(|2|p|), 



We investigate the properties of the above R.G. equations in the next section. 

IV. R.G. FLOW AND PHASE DIAGRAM 

We consider only the ferromagnetic Hund's coupling constants, so that the couplings A;,* = 1,...,5 should be 
taken to be positive in the physically relevant range. Given the initial values of Kg., — 1, . . . , 4, the R.G. flows 
are uniquely determined. The initial values are determined by Fermi momentum kp, U/t, and Ju/t- The derived 
R.G equations are valid until max[Ai/fs] ~ 0(1). If all \i converge to finite values as, t = InL oo, the initial 
fixed point is stable and the R.G. equations are valid lalong the whole R.G. trajectory since the Ai(cxD) can be made 
arbitrarily small by taking sufficiently small Jh and Im. If anij coupling constant diverges, the asymptotic behaviour 
of all coupling constants can be determined by R.G equationsQ. 

Some properties of the R.G. equations can be understood by simple inspection without numerical integration. We 
find that A2 is marginal from Eq.(p5|). The marginality of A2 is due to the chirality of interaction: it couples only 
the currents with the same chirality. The symmetric charge mode 0p+ is not renormalized as can be seen in Eq.(^0[). 
This is physically obvious since only the relative charge fluctuations affect Hund's coupling. On the contrary, we 
see from Eq.(^0|) that the Luttinger parameter of the asymmetric charge mode 0p_ renormalizes to zero when the 
R.G. equations are extended beyond its validity range. This implies that^p„ is pinned, and it is consistent with the 
formation of a charge gap from A4 renormalization (see below). In Eq.(p7l), the first term in the right hand side is 
the most dominant for the range Jh ^ U < t, in which case A4 is certainly relevant. It gives rise to both the charge 
gap of 4>p- and the spin gaps of 9a-- and 4'a-+-, as can be seen by expressing A4 term in the effective action in Abelian 
bosonized form using Eq.(^). However, for the opposite range U <C Jh < t, the first and the third terms in the right 
hand side compete, and only the numerical integrations can determine the limiting behaviour. In this range, A4 is 
expected to increase slower compared with that of the ran ge J h <^ U < t. A5 starts to grow slowly initially since the 
initial value vanishes. But eventually the first term of Eq.(^8|) dominates and A5 becomes relevant. 

The values of gaps can be estimated in the range J// <C C/ < t hy integrating the Eq.(|2^) until the A4 grows up to 
the cut-off value ta. 



t 



(31) 



t 



Note that this estimated value of the gaps is larger for the stronger Hubbard repulsion. For the range U ^ Jh < t, 
no analytical expression is available. It is worth noting that the crucial A4 interaction is not of the current-current 
interaction type, thus it is not constrained by the current conservation law and allows the presence of the anomalous 
scaling dimension. 

We have numerically integrated R.G. equations for a range of different initial values U/t, Jn/t- Some typical R.G. 
flows at C/ = and kp = 7r/4 are shown in Fig.l and Fig. 2. The Fig.l clearly shows the diverging A4 as t increases for 
a range of initial values of Jh- After the initial transient period all flows merge, which implies that the flrst term in 
Eq.(p7[) dominates. As mentioned above, for non-zero U, the divergence of A4 is much faster. In Fig. 2, the R.G flow 
for ?7 = and a very small initial value of Jh = 0.01 is shown. Initially, the A4 slightly decreases due to the last term 
in R.G Eq. (p7|) , but soon it begins to increase due to the first term of Eq.(p7|), and eventually diverges. From the 
examination of R.G fiows, we conclude that the charge and spin gaps open up for arbitrarily small initial values of Jh 
irrespective of U. Our system Eq.(^ is in the massless phase only for a line Jh = 0, which is simply two decoupled 
Hubbard chain, where each of Hubbard chain is in critical Luttinger liquid phase. In any other region with positive 
Jh > 0, A4 is relevant and the charge and spin gaps are present. 

V. DISCUSSIONS 

The spin sector of aiSingle Hubbard chain at low energy is in the same universality class of antiferromagnetic(AF) 
XXX Heisenberg chainO. Therefore, the spin sector of the system Eq.(|I]) can be approximately described as two spin 
1/2 XXX chains coupled by ferromagnetic exchange coupling at low energy. In the strong coupling limit Jh 00 
we only need to consider the triplet combination of spins from each spin chain. When projected onto this triplet 
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subspace. the coupled spin 1/2 chains becom^ssentially single spin 1 chain. The spin 1 chain is well known to be 
gappedlld, and the gap is called Haldane gapEHl. Also Strong and MilliscI studied the coupled spin 1/2 chains and 
showed that the sping gap opens even for infinitesimally small \Jh\, which is consistent with our result. The spin gap 
we have obtained can be also understood from the semiclassical point of viewEl. 

It is interesting to examine other parameter domains even though they are not directly related to our system. 
First, the antiferromagnetic exchange can be considered. It corresponds to a negative Jh, which implies the signs 
of Ai,i — 2,3,4,5 to be flipped. The R.G equation of charge Luttinger parameters, Eq.(|3^), is unchanged, but 
A3,A4,A5 now become relevant. A3 only generates a spin gap, while A4,A5 generate both the charge and spin gap. 
The relative initial values of the couplings determine the gap which is going to be opened the fastest, and R.G flows 
depend critically on the gap. The detailed investigation will not be presented here. Second, the attractive Hubbard 
interaction is another possibility (with positive Jh)- In this case, the sign of Ai is flipped and the initial value of 
is greater than 1. In this case, Ai becomes marginally relevant, and A4 becomes irrelevant. The marginally relevant 
Ai would open up a spin gap for each band I = 1,2. This is essentially the Luther-Emery spin gaplia, and the Hund's 
coupling does not play any role here. 

In this paper we have investigated two Hubbard chains coupled by Hund's coupling using renormalization group 
method for arbitrary densities away from half filling. We find that in the weak coupling regime Hund's coupling is 
always relevant, irrespective of the stength of U, and opens up gaps for both symmetric and antisymmetric spin modes 
and for antisymmetric charge mode. 
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FIG. 1. R.G. flow for three different values of Jh = 0.4, 0.2, 0.05 (Here U = and Jh is measured in units of t). 



FIG. 2. R.G. flow for Jh = 0.01 (Here U = G and Jh is measured in units of t). 
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